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ABSTRACT 

We examine the thermal and dynamical response of a neutron star to a sudden 
perturbation of the inner crust temperature. During the star's evolution, starquakes 
and other processes may deposit ^ 10'^^ ergs, causing significant internal heating and 
increased frictional coupling between the crust and the more rapidly rotating neutron 
superfluid the star is expected to contain. Through numerical simulation we study the 
propagation of the thermal wave created by the energy deposition, the induced motion 
of the interior superfluid, and the resulting spin evolution of the crust. We find that 
energy depositions of ~ 10^*^ ergs produce gradual spin-ups above the timing noise 
level, while larger energy depositions produce sudden spin jumps resembling pulsar 
glitches. For a star with a temperature in the observed range of the Vela pulsar, an 
energy deposition of ^ 10^^ ergs produces a large spin-up taking place over minutes, 
similar to the Vela "Christmas" glitch. Comparable energy deposition in a younger 
and hotter "Crab-likc" star produces a smaller spin-up taking place over ~ 1 day, 
similar to that seen during the partially time-resolved Crab glitch of 1989. 

Subject headings: stars: evolution — stars: interiors — stars: neutron — dense matter 

1. INTRODUCTION 

Glitches have now been observed in at least 20 pulsars, and are believed to be a phenomenon 
that every neutron star produces during some period in its evolution. Typical glitches in mature 
pulsars, e. g., the Vela pulsar, involve fractional jumps in the rotation rate of ~ lO^*'. Glitches in 
the relatively young Crab pulsar, on the other hand, are typically a factor of 10 to 100 smaller. 
The Crab glitch of 1989 was the first whose spin-up was partially time-resolved; following a jump 
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in the spin rate (occurring in < 2 hours), the pulsar completed the remainder of the spin-up over 
~ 1 d (Lyne, Smith. & Pritchard 1992: see Fig. 6). By contrast, the Vela pulsar has exhibited 
very different behavior in at least one case. The giant "Christmas glitch" of December 24, 
1989 (1.8 X 10"*^ fractional spin-up), which occurred during an observing session, could not be 
time-resolved; most, or possibly all, of the spin-up took place in less than two minutes (McCulloch 
et al. 1990; see Fig. 12). The Crab is unique among the known glitching pulsars in the smallness 
of its glitches. Giant glitches appear to be typical behavior among mature pulsars. No pulsar has 
exhibited any unusual behavior just prior to a glitch; these events have always occurred without 
warning. 

Several models for the origin of the glitch phenomenon have been proposed (see. e.g.. 
Ruderman 1969; Baym & Pines 1971; Pines & Shaham 1972; Pines, Shaham, & Ruderman 1972; 
Anderson & Itoh 1975; Ruderman 1976; Alpar 1977; Greenstein 1979a, b; Anderson et al. 1982; 
Alpar, Cheng, & Pines 1989; Alpar et al. 1994). In the starquake model of Ruderman (1969), 
settling of the crust under gravitational stresses abruptly decreases the star's moment of inertia, 
increasing the spin rate. At first the starquake model appeared promising, though it soon became 
clear that starc[uakes cannot accoimt for the large, frequent (about every 1000 days) glitches of the 
Vela pulsar (Baym & Pines 1971). Furthermore, since starcjuakes would occiu' over sound travel 
time scales (;$0.1 s), the starquake model is unable by itself to account for the slow 1989 spin-up 
of the Crab. A more plausible explanation is that these phenomena arise from sudden angular 
momentum transfer between the crust, whose spin rate we observe, and the more rapidly rotating 
neutron superfluirl of the interior (Anderson & Itoh 1975; Ruderman 1976). Greenstein (1979a, b) 
originally suggested that sudden heating by a starc[uake or other source could trigger increased 
frictional coupling between the crust and the core superfluid. Following the temperature increase, 
the velocity difference between the superfluid and the crust decreases; the superfluid spins down, 
while the crust spins up. However, subsequent analyses of timing noise in accretion-powered 
pulsars (Boynton & Deeter 1979; Boynton 1981; Boynton et al. 1984) indicated that ^14% of the 
neutron star interior responds to torque variations as a rigid body over time scales of 1 to 30 d, 
requiring at least part of the core to be strongly coupled to the crust. In support of this constraint, 
theoretical study by Alpar, Langer, & Sauls (1984) suggests that the entire core is tightly coupled 
to the charged components of the star. In light of these developments, glitches are now believed 
to originate in the inner crust, where interactions between superfluid vorticity and nuclei create a 
metastable state in which the superfluid is weakly coupled to the more slowly rotating crust. 

Here we consider the problem of the thermal and dynamical response of a neutron star to a 
sudden local increase of the inner crust temperature. We employ the results of recent studies of 
the frictional coupling between the inner crust superfluid and the crust, which indicate that the 
coupling is a strongly increasing function of temperatiu'e (see. e. g.. Alpar 1977: Alpar et al. 1984; 
Alpar, Cheng & Pines 1989; Link & Epstein 1991, hereafter LE; Link, Epstein, & Baym 1993, 
hereafter LEB; Chau & Cheng 1993a, b). We find that energy depositions of ~ 10^'^ ergs produce 
gradual spin-ups above the timing noise level, while larger energy depositions produce sudden spin 
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jumps resembling pulsar glitches. For a star with a temperature in the observed range of the Vela 
pulsar, an energy deposition of ~ 10^^ ergs produces a large spin-up taking place over minutes, 
similar to the Vela "Christmas" glitch. Comparable energy deposition in a younger and hotter 
"Crab-like" star produces a smaller spin-up taking place over ^ 1 day, similar to that seen during 
the partially time-resolved Crab glitch of 1989. 

The paper is organized as follows. In §2 we summarize the current understanding of the 
neutron star interior and give an overview of the problem we are considering. In §3 we discuss 
the physical conditions of the neutron star inner crust. In §4 we present a descriptions of the 
superfluid, crust and thermal dynamics. In §5 we describe our numerical simulations of the spin-up 
process. In §6 we discuss our results, and in §7 we present our conclusions. 

2. OVERVIEW 

2.1. The Physical Setting 

A neutron star consists of ~ 1.4Mq of mostly neutrons in /3-equilibrium. Beneath a solid 
or liquid surface, begins the outer crust which contains a lattice of nuclei and relativistically 
degenerate electrons. At a density of --^^ 4 X 10"'^"'^ g cm^"^ the inner crust begins, consisting of an ^Sq 
neutron superfluid, a lattice of neutron- rich nuclei and a sea of relativistically degenerate electrons. 
Near the density of nuclear matter, /9„„, = 2.8 x 10^^ g cm^'^. the inner crust dissolves into a 
fluid core consisting of mainly "^^2 superfluid neutrons and a small admixture of superconducting 
protons with normal electrons. Depending on the the behavior of matter above nuclear density, 
and the mass of the star, there may exist a distinct inner core containing condensed pions, quarks 
or some other exotic state. 

The neutron superfluid is expected to play a crucial role in the star's rotational and thermal 
evolution. As a magnetized spinning neutron star slowly loses angular momentum, the external 
torque acts directly on the charged components, which include the crustal lattice, electrons and 
protons. The inner crust superfluid, however, couples only weakly to the charged components 
and maintains a higher rotation rate than that of the crust. Angular momentum transfer from 
the superfluid to the crust, whether sudden or gradual, affects the star's spin dynamics. Heat 
generated by the differential rotation between the superfluid and normal components affects the 
star's thermal evolution (Greenstein 1975: Harding, Guyer. & Greenstein 1978: Alpar et al. 1987; 
Shibazaki k Lamb 1989; Van Riper 1991; Van Riper. Epstein, k Miller 1991: Umeda et al. 1993; 
Van Riper, Link, k Epstein 1995). Coupling between the long-term thermal evolution and the 
rotational dynamics has been considered by Sliibazaki k Lamb (1994) and Sliibazaki k Mocliizuki 
(1994). 

The rotating superfluid is threaded by an array of vortex lines, aligned with the star's rotation 
axis. The vortex arrangement determines the rotational state of the superfluid. Within each 
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vortex line core, the fluid is normal. Processes by which a vortex interacts with the charged 
components include the scattering of electrons, phonons, impimties and lattice defects with 
thermally excited neutrons in the vortex core (Feibelman 1971: Harding. Guyer, & Greenstein 
1978; Jones 1990b), scattering of electrons with the charge distribution induced about the vortex 
core by its interactions with nuclei (Bildsten & Epstein 1989), coupling of vortex translational 
motion with the electron-phonon system of the lattice (Jones 1990a), and the coupling of vortex 
Kelvin modes (kelvons) with lattice phonons (Epstein & Baym 1992: Jones 1992). For angular 
momentum exchange between the superfluid and the crust to occiu' through any of these processes, 
the vortices must change their arrangement. In the inner crust, however, it is expected that 
interactions between the ^Sq superfluid vortices and the lattice nuclei pin the vortices to the 
lattice, nearly flxing the superfluid angular velocity and allowing a significant velocity lag to 
develop between the superfluid and the crust as the star slows down (Anderson & Itoh 1975; 
Alpar 1977; Alpar et al. 1984; Epstein & Baym 1988). The superfluid flowing past a pinned vortex 
lines exerts a radially outward force on it. In the presence of thermal fluctuations, a pinned 
vortex line occasionally unpins and creeps outward (see, e.g., Alpar et al. 1984; LE; LEB; Chau 
& Cheng 1993a, b), the superfluid loses angular momentum and exerts torque on the crust. The 
coupling between the superfluid and the crust through thermally-activated vortex creep is highly 
temperature dependent. If the coupling suddenly increases, from a jump in the temperature, for 
instance, angular momentum transfer to the crust could drive a spin-up.^ 

In contrast to the inner crust superfluid, the core superfluid may be tightly coupled to the 
charged components. The vortices of the core are strongly magnetized due to entrainment of the 
superfluid protons by the neutron circulation associated with each vortex (Alpar, Langer & Sauls 
1984). Electrons efficiently scatter off the vortex line magnetic moments, coupling the charged 
components and superfluid together, effectively creating a single fluid. The fluid core is likely 
coupled to the crust by magnetic stresses (Abney, Epstein, & Olinto 1995). 



2.2. The Thermally-driven Spin-up 

We conjecture that glitches originate from the sudden deposition of thermal energy in the 
inner crust. Because the frictional coupling between the crust and the inner crust superfluid 
is highly temperature dependent, any process that leads to inner crust heating will cause a 
substantial increase in the coupling. As the star evolves, internal heating can occur through 
a number of processes, including, internal friction, accretion, plastic or sudden relaxation of 
structure, and nuclear reactions. Starquakes are expected to occur as the star slows, the centrifugal 
force at the surface diminishes, and the rigid crust settles. In addition, as the star slows, magnetic 



^ Jones (1990c; 1991) has considered the possibility that vortex lines do not pin, and that the superfluid maintains 
a state of near corotation with the crust. In such a case the superfluid could never acquire a sufficient velocity 
difi^erence from that of the crust to drive glitches. 
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flux may be swept from the core by the outward migration of vortex lines. Even though the 
classical ohmic dissipation time scale is > 10^*^ years for magnetic flelds threading the core of a 
neutron star (see, e. g., Urpin & Van Eiper 1993), flux tubes may be expelled much more rapidly. 
Magnetic pressure in the flux tubes makes them lighter, causing them to float upward (see, e. g., 
Jones 1987). Another possibility is that the flux tubes are dragged out by the radial motion of 
neutron-superfluid vortex lines as the star slows. (Srinivasan et al. 1990; Ding, Cheng & Chau 
1993). In either case the magnetic field topology may evolve in such a way as to produce an abrupt 
change in direction and strength near the core-crust interface. The magnetic stresses resulting in 
the crust could induce quakes. 

For deflniteness. we envision the heating occurring suddenly during a starquake. We estimate 
the average accumulation rate of the strain energy from gravitational stresses as (Baym & Pines 
1971; Cheng et al. 1992) 

BO,t-^„ (1) 



dt 

where B (~ 10^*^ ergs) is the crust strain energy. 0^. is the critical strain angle at which the lattice 
breaks, and ^age is the pulsar age. The value of 6c is uncertain, and estimates range between 10~^ 
(characteristic of the Earth's crust) to 0.1 for a perfect Coulomb lattice. If the accumulated strain 
energy is released every few years, the Crab and Vela could experience energy depositions in excess 
of 1042(61^/10-2) ergs. 

Following the energy deposition, a thermal wave propagates through the inner crust, the 
superfluid of the ever-growing heated region transfers angular momentum to the crust, and the 
crust spins up. The glitch ends when the first of two things happens: either the thermal energy 
becomes highly diluted over a large volume, or the superfluid of the active region nearly achieves 
corotation with the crust. The magnitude and time scale of the spin-up are determined by the 
relative temperature change caused by the energy deposition. The larger the heating, the greater 
the coupling rate and the faster the spin-up. Also, the greater the energy deposition, the larger 
the region that is ultimately affected, and the larger the resulting spin-up. 

In the context of this physical picture, we make two remarks: 

• Spin jumps begin suddenly. The sudden heating that accompanies a starquake causes the spin-up 
to begin without warning. 

• For a given energy deposition, cooler pulsars produce larger, fa.ster glitches. As a neutron star 
cools, its specific heat drops. In a cooler star, with its relatively low specific heat, a given energy 
deposition gives rise to a relatively large change in the temperature, giving rise to a greater 
frictional coupling than in the case of a hotter pulsar; the glitch is correspondingly larger and 
faster. 

The first statement is consistent with all glitches detected so far. The second remark is 
suggestive of the cause behind the markedly different glitch magnitudes and time scales seen 
during the 1989 Crab glitch and during the Vela Christmas glitch. Before turning to a quantitative 
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description of the thermally-driven spin-up process, we discuss the physical conditions that prevail 
in the inner crust. 



The nuclear matter equation of state (EOS) determines the bulk properties of a non-rotating 
neutron star, including the radius, total moment of inertia, inner crust thickness, and inner crust 
moment of inertia. Recent calculations based on the best potentials currently available for the 
nucleon-nucleon interaction and including three-body interactions give EOSs that are moderately 
stiff to stiff (see, e. g., Wiringa. Fiks. & Fabrocini 1989). Moreover, determination of the minimum 
mass of the inner crust supcrfluid needed to account for post-glitch relaxation behavior (Link, 
Epstein, & Van Riper 1992), considered along with the mass measurement of PSR 1913-1-16 
(Wolszczan 1991), rules out very soft EOSs. For the purposes of determining the effective extent 
and characteristic density of the inner crust superfluid, we will consider the model of Friedman & 
Pandharipande (1981; FP), a representative EOS of moderate stiffness. 

Fig. 1 shows the inner crust density as a fimction of radius for a IAMq neutron star built 
on the FP EOSs. The inner crust ranges in density from 4 x 10^^ g cm"*^ to near nuclear matter 
density. /9„„, = 2.8 x lO' ^ g cm^'\ Recent calculations indicate that the the lattice dissolves at a 
density ^ O.C/9„to (Lorenz, Ravenhall, &: Pethick 1993; Pethick, Ravenhall, & Lorenz 1995). Most 
of the moment of inertia of the crust is in the region of highest density, and it is in this region that 
dynamic coupling between the pinned superfluid and the crust can give the largest contribution 
to variations in the spin rate of the crust. About 90% of the inner crust mass is in the region 
9.8 < r < 10.0 km with a mean density of p 1.5 x 10"^^ g cm"*^. We take this average density as 
characteristic of the inner crust superfluid. 



Soft x-ray emission has been detected from at least four neutron stars. If Tg is the star's 
measured effective surface temperature, we can estimate the internal temperature, T, using the 
following expression appropriate for a non-magnetized star (Gudmundsson, Epstein, & Pethick 



3. 



PHYSICAL CONDITIONS IN THE INNER CRUST 



3.1. The Equation of State 



3.2. Temperature 



1982) 




(2) 



where 
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and 

e-^=(l-^)"''. (4) 

Here Tg = T^g = T^/IO*^, gsi4 = (/^/lO' ' is the surface gravity, M and R are the stellar 

mass and radius, and Rgch is the Schwarzchild radius. 

The effective surface temperature is related to the luminosity seen by a distant observer 
through, 

AirR^aT^ = e-'^^Loc. (5) 

In some cases where surface emission has not been detected {e.g., the Crab pulsar), upper limits 
on Loo are available, and we can obtain an upper limit on the internal temperature using eqs. 
[2]-[5]. For the purpose of estimating internal temperatures, we identify Tg with the temperature 
obtained from a black body fit to the observed portion of the spectrum. The dependence of the 
atmospheric opacity on photon energy, however, can lead to a difference between the temperature 
inferred from such a blackbody fit and the true effective temperature (see, e.g., Romani 1987; 
Miller 1992, 1993; Shibanov et al. 1993). This issue has yet to be fully resolved. 

Fig. 2 summarizes our knowledge of neutron star internal temperatures. Four detections and 
one upper limit exist in the age range 10^ — 10^ years. These data suggest that stars in this age 
group have similar internal temperatures. 



4. DYNAMICAL DESCRIPTION 

We consider a star of total moment of inertia /, slowing under the influence of a constant 
external torque, IO.oc- The total angular momentum J is given by the angular momentum of the 
crust and all components strongly coupled to it, plus the angular momentum of the superfluid Jg, 

J(t) = Jc{t) + Js{t) = iMt) + J d^r r'^pMr,t) = Jq - I|fioo|i, (6) 

where Jq is the total angular momentum at t = 0. We take to be the moment of inertia of the 
solid crust plus that of the fluid core. The star is in a rotational steady state if the crust and 
superfluid are both slowing at the same rate, Qoo- 



4.1. The Superfluid 

A rotating superfluid is threaded by an array of vortex lines aligned with the rotation axis, 
each contributing a quantum of vorticity k to the net superfluid flow. [For a neutron superfluid, 
K = /i/2m„ where m„ is the neutron mass]. For a superfluid with an average flow velocity Vgi''') 
and angular velocity Qg(r), the total circulation is given by the total number of vortex lines within 
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an integration contour enclosing the axis of rotation, 



Vg ■ dl = 2TTr'^Qs{r) = k. dr'2'Kr'n{r') 

Jo 



(7) 



where r is measured from the rotation axis and n{r) is the local areal density of vortex lines. 
Differentiating c.q. [7] with respect to r, we obtain a differential equation relating the vortex 
density to the local superfluid angular velocity, 



Radial vortex motion changes the number of vortex lines within a radius r. and hence the 
superfluid rotation rate (see c.q. 7). For vortex lines crossing a cylindrical boundary at r with a 
radial velocity component w(t), number conservation requires 



where N is the total number of vortex lines contained within the volume. An equivalent expression 



Eq. [10] is the basic dynamical equation describing the superfluid rotation. The radial component 
of the vortex velocity, v. is determined by the manner in which pinned vortices impin and move 
through the inhomogeneous environment of the crustal lattice, a problem which has been studied 
in various limits (LE; LEB; Chau & Cheng 1993a, b). We use here the results of LE and LEB, 
which have the most general applicability. A pinned vortex segment unpins by overcoming an 
activation barrier whose height decreases as the lag velocity between the superfluid and the crust, 
w = Jl, — Qc, is increased. The mechanics of the impinning process, and hence, the height of the 
activation barrier, is determined by u) and by the importance of vortex tension relative to the 
pinning force per nucleus. Estimates of the vortex stiffness (LEB) based on the gap calculations 
of Ainsworth et al. (1989) indicate that vortices are very stiff in the densest regions of the inner 
crust, such that unpinning of a vortex segment in these regions requires the breaking of multiple 
pinning bonds simulaneously. As the lag u is reduced, more pinning bonds must break, and the 
activation energy grows as cj"^. For the neutron star internal temperatures of interest here, stiff 
vortices overcome their pinning barriers primarily through classical thermal activation, rather 
than through quantum tunneling. We use the following expression for the thermally-activated 
vortex creep velocity (eq. 6.9, LEB with eq. B.12, LE): 




(8) 




(9) 



IS 




(10) 



v{u), T) = vq e 



13/Tu, 



(11) 



This form applies in the limit w -C cjc, where ujc is the critical lag at which vortices cannot remain 
pinned: our use of this limit is justified below. The quantity /3 measures the strength of the 
coupling between the superfluid and the normal matter; it depends only on the local mass density. 



- 9 - 



and is the same for a region of a given density for all stars. The coupling parameter can be 
written in terms of and the effective pinning strength Ep per nucleus, and can also be related 
to observable quantities (eqs. 6.11 and 6.14, LEB): 

13 = 0.54 uj.Ep \nMn)n\n^ J^4«o^^ ^ ^^2) 

where ^age = f^s/2|r2oo|- Here Tq is the temperature in the absence of sudden heating and tr is 
the coupling time between the superfluid and the normal matter. The dependence of p on fage 
and Tq in eq. [12] is only apparent: t,, scales as (|Ooo|Tb ln^[4wotage/-^]) "'^ (LEB). For a relaxation 
time measured from post-glitch spin behavior and an internal temperature deduced from surface 
thermal emission, p can be determined. The exact value of vq is unimportant, and is expected to 
be ~ 10"^ cm s^^ (LE; Epstein & Baym 1992; LEB). 

If the inner crust temperature is perturbed above Tq, the vortex velocity increases dramatically, 
the superfluid rotation rate decreases (see eq. 10), and the crust spins up. As the superfluid 
approaches corotation with the crust, the vortex velocity drops precipitously, and the superfluid 
is essentially decoupled horn the crust. Immediately after the glitch, the external torque acts on a 
lesser moment of inertia, and the crust spins down at a greater rate. 



4.2. Thermal Dynamics 

The time-dependent coupling of the superfluid to the crust is determined by the evolution of 
the local temperature. Following the sudden energy deposition, a thermal wave moves through the 
crust according to 

C.^ = V-(ktVT), (13) 

where and kt are the local specific heat and thermal conductivity, both functions of 
temperature and density. In eq. [13] we neglect the heat generated by the friction between the 
two components; as discussed in §6, heat generation is unimportant for stars in the temperature 
range of the Crab and Vela, but is likely to be important for colder stars. The specific heat resides 
mostly in the degenerate electrons and scales linearly with temperature. Thermal conduction is 
inhibited by phonons, and the conductivity decreases with increasing temperature. 



5. SIMULATIONS: GEOMETRY AND INITIAL CONDITIONS 

We model the portion of the inner crust superfluid that drives the spin-up as a cylindrical 
shell. The greatest contribution to the moment of inertia of the inner crust superfluid is in a 
band about the equator of the star, and a cylindrical geometry simplifies our calculations while 
capturing the essence of the inner crust dynamics. Furthermore, we expect crust cracking in a 
slowing neutron star to occur mainly near the equator, where the gravitational stresses become 
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the greatest. A full reproduction of the many aspects of pulsar timing behavior is beyond the 
scope of our simple model. The multi-component nature of post-glitch recovery is presumably due 
to the contributions of distinct inner crust regions which find themselves in different rotational 
states following the glitch (see, e. g., Alpar et al. 1984). We will focus here instead on reproducing 
the spin-ups themselves, which occur more rapidly than the post-glitch recovery. 

The effective crust in our model is a cylindrical shell of constant supcrfluid mass density 
p, thickness A, and height /i, with outer and inner radii i^c ± A/2 (see Fig. 3). We define a 
cylindrical coordinate a;, such that a; = at radius R^. We estimate the height of the effective 
pinning region by the height at which at a vortex a distance i?,^. from the axis exits a spherical 
shell of radius i?,^. -|- A/2. For an FP EOS. ft, is ~ 0.31?,^-. The height of the heated region is more 
difficult to estimate, but could be comparable to that of the pinning region, as we assume. 

We consider the deposition of energy E in the cylindrical effective crust, such that the greatest 
heating is at the center of the cylinder, and the temperatiu'e falls off towards its inner and outer 
boundaries. The heating is independent of position along the rotation axis. The temperature 
distribution is in a gaussian profile of half- width a in the x direction, «. e., 

AT(a;,i = 0) = AT(0,Q)e-^'/2'"'. (14) 

For a specific heat that depends linearly on T, the energy deposited is given by 

/i/2 fT 



E = 2Tr dz dr dT C„(T) 

J-h/2 Jr^-A/2 JTn 



l-h/2 JRc-A/2 JTo 

f°° , (AT(x.t = (]) 1 /AT{x.t = (])\-] , , 

^ 2nR,hCM dx I '- + - i , (15) 

where the subscript denotes quantities evaluated at the unperturbed temperature, Tq. We 
take cr -C A, and have extended the integration over an infinite domain as a simplification. 
Additionally, we neglect the effects of thermal diffusion into regions of lower density. 

The solution to eq. [15] is 

'"^<»'"'- .v^+('2 + -,,^-? (16) 



To V ■K^/'^CMRcha 

We are interested in situations where a <^ h < R, so considering one-dimensional diffusion of a 
planar wave in the x direction is a good approximation. In this preliminary study we neglect 
the effects of the temperature dependence of the thermal conductivity and specific heat on the 
evolution of the diffusion wave, and solve the diffusion eq. [13], which is now 

dt Cm ■ ^ ' 

This linear approximation is rigorously valid only for temperature enhancements significantly 
less than unity; however, below we show this approximation to be adequate for temperature 
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enhancements of order unity as well. The solution to eq. [17] with our initial conditions is 

AT{x,t) = AT(0,0) (^j^ ^^\-toxV2ta^^ ^^g^ 

where tQ = C^qg'^ I2kq. 

Assuming that the superfiuid and the crust are in rotational eciuilibrium before the energy 
deposition, both components slow at the same rate iloo- To obtain the initial lag throughout the 
crust, we solve eqs. [10] and [11] numerically with Qs{r,t) = —\Q.oc\- We then deposit an amount 
of energy E, and solve eqs. [16] and [18] for the temperature throughout the crust as a function of 
time. At each time step, we solve eqs. [10] and [11] for Cls{r,t), eq. [6] for Clc{t), and then update 
u{r,t). 



6. DISCUSSION 

6.1. Example Simulations 

To illustrate the consequences of sudden inner crust heating, we present solutions for three 
models (see Tables 1 and 2). For the first model, the parameters were chosen as representative of 
the Crab pulsar. The second model represents Vela, and the third describes an older star, 5 x 10^ 
years in age. The Crab's surface temperature is taken to be about half the observational upper 
limit. The temperatiu'es of Vela and the older star are the same, about one-third that of the Crab. 
The superfluid coupling parameter depends on the strength of vortex pinning only, and to the 
extent that spin-up events in different pulsars are driven by regions of the same density, we expect 
/3 to be a constant among pulsars; we take /3 = 2 X 10^ rad s^^ K in each of our three models. For 
this choice, the superfluid relaxation times implied by eq. [12] for our assumed temperatures are 
495 d for Vela and 7 d for the Crab, consistent with the post-glitch relaxation times seen in 
these two pulsars. The lag velocities in the steady state are 0.04 rad s^^ for the simulation of the 
Crab, and ^0.12 rad s^^ for the simulations of Vela and an older pulsar. As these lag velocities 
are signicantly less than the estimates of Uc — 1 rad in the densest regions of the crust (LEB, 
Table 1, using the gap parameters of Ainsworth et al. 1989), our use of the low-lag form for the 
vortex velocity (eq. 11) is justified. The stellar radius, moment of inertia, effective crust extent, 
and average superfluid density are based on the FP EOS. For an FP star, 90% of the crust's mass 
is contained in the region 9.8 < r < 10.0 km; we take A, the thickness of the effective crust to be 
200 m, and R^, the radius of the center of the effective crust, to be 9.9 km. The initial width of 
the heated region is cr = 20 m for each model. We now discuss example simulations from each of 
our three models. 

Crab- like Simulation. As our first example, we consider an energy deposition of 2.1 X 10^^ 
ergs. The progress of the thermal wave during the glitch is shown in Fig. 4. The energy 
deposition produces a two-fold increase of the initial inner crust temperature. As the thermal wave 
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propagates, the superfluid in the heated region loses angular momentum (see Fig. 5). The angular 
momentum lost from this region is given to the crust, whose spin behavior we show in Fig. C. 
About half of the angular momentum transfer occurs within the first 100 seconds. The remainder 
of the spin-up occurs over the following day, and amounts to a final spin-up of AQc/Qc — 7 x 10~^. 
After ~ 1 d, the thermal energy has been diluted through the crust, and the superfluid relaxes 
over the following days, during which the crust spins down at a greater rate than before the glitch. 
Dimng this time, the superfluid in the region surrounding x = has reduced the lag to the degree 
that vortex motion is effectively turned off. the superfluid of this region is decoupled from the 
crust, and the external torque acts on a lesser total moment of inertia. The temperature change 
produces only small changes in the local vortex density, depletions and accumulations of less than 
1% (see Fig. 7). In spin-up magnitude and time scales, as well as post-glitch relaxation, this 
simulated spin-up is similar to the 1989 Crab glitch. 

Vela-like Simulation. This model differs from the model for the Crab in the unperturbed 
temperature (about three times lower), and the star's angular velocity and spin-down rate. For 
this example, the energy deposition is 1.51 X 10^^ ergs. Because of the lower specific heat of this 
cooler pulsar, the energy deposition produces a nearly six-fold increase in the initial temperature 
(Fig. 8), resulting in a significantly larger and faster change in the superfluid rotation rate 
(Fig. 9). Roughly half of the angular momentum transfer occurs within the first second, with 
the remainder occurring over the subsequent several hundred seconds, giving a final spin-up of 
AQc/^c — 2 X 10~^ (Fig- 10). Significant vortex motion occurs, with vortex depletions and 
accumulations of order unity (see Fig. 11). This simulated spin- up is consistent in both magnitude 
and time scale with the Vela Christmas glitch. In Fig. 12 we show the predicted pulse arrival 
times, compared to the data from the Christmas glitch. Lesser energy input produces a smaller, 
slower glitch; an example is shown in Fig. 13. Small events of this short have been seen in Vela. 

An Older Pulsar. In Fig. 14 we show the response of an older, more slowly decelerating 
star to an energy deposition of 9 x 10^^ ergs. For an energy deposition of 1.51 x 10^^ ergs, as 
considered in the Vela-like simulation, the spin-up is ~ 2.5 times larger, over a shorter time scale. 

An interesting question is how the glitch magnitude and time scale depend on the energy 
deposition; these quantities are shown in Figs. 15-17. Spin-ups above the timing noise level 
(AQc/f^c ~ 10^^^) occur for energy depositions above ~ 10^*^ ergs. For energy depositions ~ 10^^ 
ergs and larger, a given energy deposition produces a spin-up nearly two orders of magnitude 
larger in the Vela-like simulation than in the Crab-like simulation. In the simulation of an older 
pulsar, the behavior at low energies is similar to that of the Crab. For large energy deposition, 
however, the resulting spin-up is somewhat larger than in the model of Vela. For very large energy 
depositions (~ 10^*^ ergs) the curves in Fig. 15 begin to flatten out; in these cases the energy 
deposition causes nearly the entire inner crust superfluid to achieve corotation with the crust, 
giving the largest possible glitch. Energy depositions of such magnitude are signiflcantly greater 
than those expected from starquakes. In Figs. 16 and 17 we show the fraction of the total spin-up 
after selected times as a function of energy deposition. The rate of the spin-up increases markedly 
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with energy, and hence, with ghtch magnitude. 



6.2. General Properties of the Solutions 

Dependence on temperature. The quantity that most strongly affects the spin-up 
magnitudes and time scales among our simulations is ChqTq i"x Tq . which determines the relative 
temperature change for a given energy deposition (see o.q. 16): it is a factor of ^ 10 lower in 
our Vela-like and old pulsar simulations than in the Crab-like simulation. Consequently, a given 
energy deposition produces a larger (up to a factor of 10) temperature change in the Vela-like and 
older pulsar simulations, leading to a larger and faster transfer of angular momentum from the 
superfluid to the crust. The glitch is correspondingly larger and faster. 

Dependence on the superfluid coupHng parameter. The glitch magnitude depends 
sensitively on the magnitude of the superfluid coupling parameter (3. For a given value of /3 the 
lag in the steady state, prior to the energy deposition, is from eqs. [10] and [11] 

.„.|.,„-il!5|ss£. (19) 

IQ K 

For larger /3, the superfluid is less tightly coupled to the crust and ojss is larger. After the energy 
deposition, the local vortex velocity becomes 



v{u), T) = vq exp 



(20) 



where Suj is the local change in the lag from ujss as the spin-up proceeds (Suj < 0). For larger /3 and 
uigs, the vortex velocity decreases less rapidly as the lag decreases, and the glitch is correspondingly 
larger. 

Dependence on pulsar age. An older pulsar has a lower vortex creep velocity in the steady 
state and a lower uigg. The vortex velocity in an older star increases less with the temperature 
increase, and drops more rapidly as the lag decreases (see cq. 20). For small or modest energy 
input, the spin-up is smaller (compare the Vela-like and old pulsar simulations in Fig. 15). For a 
pulsar of any age, however, a very large energy deposition brings the entire inner crust superfluid 
to a state of near corotation with the crust, giving rise to a glitch of magnitude 



(21) 



where Ig is the inner crust moment of inertia. The glitch magnitude in an old pulsar is now larger 
than in a younger pulsar of the same temperature by a factor of up to [The decrease in 

LOss with pulsar age is relatively small]. In our model of an old pulsar, the spin rate is ~ 7 times 
smaller than in our model of Vela, corresponding to a glitch larger by the same factor for energy 
inputs of ~ 10^*^ ergs (see Fig. 15). 
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Dependence on dimensions of initial heating. Tlio total spin-up and the time scale over 
which it occurs are rather insensitive to the initial width of the energy deposition, a; the glitch 
magnitude changes by a factor of 2 as c is varied from 5 m to 100 m. The glitch magnitude 
increases slowly with the height of the effective region, h. 

6.3. Other Effects 

Non-linear diffusion. The large 6-fold) increases in the local temperature we have 
considered in our Vela-like simulation raise the question of the adequacy of our linear treatment 
of the diffusion problem. Accounting for the temperature dependence of the heat transport would 
quantitatively change our results. To estimate the importance of these non-linear effects, we have 
performed a simulation with our Vela-like model using a thermal diffusivity that is smaller by a 
factor of ~ 8, corresponding to a temperature ~ 3 times higher. The initial 300 seconds of the 
evolution closely follow the results of the original simulation at lower temperature. Owing to the 
lower thermal diffusivity, the final portion of the spin-up occurs over ~ 1 hour. The resulting 
spin-up is about 25% larger because more vortices can move before this more slowly moving 
thermal wave decays. We expect the results presented here to be qualitatively unchanged by the 
real effects of non-linear diffusion. 

Frictional heating. In this study we have ignored the back-reaction of the frictional heating 
that accompanies the glitch. This heating further increases the frictional coupling. An important 
question to be resolved is if a feedback process could lead to a thermal runaway, as suggested 
by Greenstein (1979a, b), wherein an enormous spin jump occurs, accompanied by a "thermal 
explosion" from the stellar surface. The effects of glitch-induced heating become important when 
the heat generated is comparable to the initial energy deposition. For a lag uj prior to the glitch, 
the heat liberated is A.Tw. where A, 7 is the change in angular momentum. The lag depends on 
the strength with which vortices pin to the inner crust lattice, and we estimate that it is ^ 10 
rad s~^ (LE)- In the Vela pulsar A J is of the order of 10^^ erg rad~^ s, typically, and the heat 
generated during the glitch is ^ 10^^ ergs, less than the energy deposition considered in our 
example simulation. However, we expect that in a cold pulsar, with its correspondingly lower 
specific heat, giant glitches can arise from energy depositions significantly less than 10^^ ergs; we 
anticipate the effects of frictional heating to play an important role in the rotational and thermal 
dynamics of these pulsars. We will consider this possibility in future work. 

7. CONCLUSIONS 

As a neutron star evolves, any process that leads to internal heating will affect both the 
thermal and rotational evolution of the star. In particular, starc^uakes could produce significant 
internal heating, dramatically increasing the frictional coupling between the crust and the more 
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rapidly rotating neutron superfluid of the interior. In this paper we have studied the thermal and 
dynamical response of a neutron star to a sudden deposition of energy in the inner crust, as would 
result from a starquake. We find that energy depositions of ~ lO '*^ ergs produce gradual spin-ups 
above the timing noise level, while larger energy depositions produce sudden spin jumps resembling 
pulsar glitches. The differences in spin behavior between the Crab pulsar and other, more mature 
pulsars may be attributable to differences in temperature. For a star with a temperature in the 
observed range of the Vela pulsar, an energy deposition of ~ 10*^^ ergs produces a large spin-up 
taking place over minutes, similar to the Vela "Christmas" glitch. Comparable energy deposition 
in a younger and hotter "Crab-like" pulsar produces a smaller spin-up taking place over ~ 1 day, 
similar to that seen during the partially time-resolved Crab glitch of 1989. 
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Fig. 1 — Inner crust density versus distance from the stellar center for a star constructed from 
an FP EOS. The dashed line indicates a density of Q.Gpnm, at which the calculations of Lorenz, 
Ravenhall, & Pethick (1993) indicate the nuclei dissolve. 

Fig. 2 — Neutron star internal temperatures. Note the similarity in temperatures in the age range 
Iq4 _ iq6 

years. Except for the Crab pulsar, for which the age is known, an n = 3 breaking index 
model was used to obtain the pulsar age. An FP EOS was assumed. 

Fig. 3 — Model geometry. The effective crust is modeled as a cylindrical shell. Heating occurs in a 
region of height h and thickness 2a. 

Fig. 4 — Evolution of the inner crust temperature enhancement in our Crab-like simulation, 
following the deposition of 2.1 X lO '" ergs. 

Fig. 5 — Evolution of the superfiuid. We show the change in the superfluid angular velocity, 
relative to its value before the energy deposition. 

Fig. 6 — The response of the crust following the deposition of 2.1 x lO '- ergs at t = 0. The 
pre-glitch spin evolution has been subtracted. The data from the 1989 glitch are shown for 
comparison. About half of the spin-up occurs during the first two minutes. At late times the 
superfluid is decoupled from the crust, the external torque acts on a lesser moment of inertia, and 
the spin excess decays. 

Fig. 7 — Evolution of the vortex density, in units of the initial density no 

Fig. 8 — Evolution of the inner crust temperature enhancement in our Vela-like simulation, 
following the deposition of 1.5 x 10^^ ergs. 

Fig. 9 — Evolution of the superfluid. 

Fig. 10 — The response of the crust following the deposition of 1.5 X 10^^ ergs at f = 0. The 
pre-glitch spin evolution has been subtracted. About half of the spin-up occurs (hiring the first 
three seconds. Decay of the spin excess, not visible on this scale, occurs at late times. 

Fig. 11 — Evolution of the vortex density. 

Fig. 12 — The pulse arrival times for the simulation of Fig. 10. The arrival times expected from 
the pre-glitch behavior have been subtracted. Data from the Christmas glitch are shown. 

Fig. 13 — A small, slow spin-up and the subsequent recovery resulting from the deposition of 
2 X 10^'^ ergs in our model of the Vela pulsar. 

Fig. 14 — The response of an older star (age 5 x 10^ years) to the deposition of 9 x 10*^^ ergs at 
t = 0. The unperturbed temperature is equal to that of our Vela-like simulation. 

Fig. 15 — Total spin-up as a function of energy deposition. 

Fig. 16 — Fraction of the total spin-up versus glitch magnitude at selected times, for our model 
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of tho Crab pulsar. The upper and lower data points correspond to the total spin-ups observed 
during the 1989 Crab glitch after 1 d and 0.1 d, divided by the spin-up achieved after 1.64 d (the 
highest data point appearing in Fig. 6). 

Fig. 17 — Fraction of the total spin-up versus glitch magnitude at selected times, for our model of 
the Vela pulsar. 



- 18 - 



FP Inner Crust 




r (km) 



Figure 1 



- 19 - 



CD 



10 



C/) 

b 

a5 



10 







10 



-1 



Crab 



V 



Vela 0656+14 



1916+14 



1 055-52 



Geminga 



1929+10 - 



Y 



10' 



10^ 10^ 10^ 10^ 10^ 



Pulsar age (years) 



Figure 2 



- 21 - 



Crab-like Simulation 



t=1 s 

t=10's 




x(m) 



Figure 4 



- 22 - 



Crab-like Simulation 




100 



x(m) 

Figure 5 



Crab-like Simulation 




Figure 6 



- 24 - 



Crab-like Simulation 



1.010 



1.005 



1.000 - 



0) 
"D 

X 
CD 



O 

> 0.995 




0.990 



-100 



■50 





x(m) 



50 



100 



Figure 7 



- 25 - 



Vela-like Simulation 




x(m) 



Figure 8 



- 26 - 



Vela-like Simulation 



- 



-2 - 



o 

CO 



CO 



a 



-4 - 



1- -6 - 



-8 - 




\ \ \ 



\ \ 
\ ' I 



W 

w 



W 



/ 

/ / 
/ / / 



/ // 

/ // 
/ '/ 

/ I' 
' // 

/ ; 

a - 

'/ 



t=0.01 s J 
t=10s 
t=10's 
t=2x1 0' s 



11 
I 



-10 



-100 



■50 





x(m) 



50 



100 



Figure 9 



- 27 - 



Vela-like Simulation 



200 



-| ' 1 ' r 



1 ' 1 ' r 



150 - 



o 100 

a 

a 

< 



- 



■50 



_i \ 1 L 



J 1 L 



J 1 \ u 



■20 -10 10 20 30 40 50 60 70 



time (minutes) 



Figure 10 



- 28 - 



Vela-like Simulation 



1.3 



1.2 - 



t=0.01 s 

t=10s 

t=10's 

t=2x10's 



1.1 - 



CD 
"D 

X 
0) 

o 
> 



1.0 - — 



0.9 - 



0.8 - 




0.7 



-100 



■50 





x(m) 



50 



100 



Figure 11 



Vela-like Simulation 




Figure 12 



- 30 - 



Vela-like Simulation 



0.15 



0.10 



o 



a 



o 



a 



0.05 - 



00 



0.00 - 



-0.05 




-500 



500 1000 1500 2000 



time (days) 



Figure 13 



- 31 - 



Older Star 



200 



150 - 



o 100 

a 

a 

< 

50 



- 



■50 




-5 



5 10 15 

time (minutes) 



20 



Figure 14 



- 32 - 




Energy deposition (ergs) 



Figure 15 



- 33 - 



Crab-like simulation 




Figure 16 



- 34 - 



Vela-like simulation 




c c 



Figure 17 



- 35 - 



TABLE 1 

MODEL SIMULATION PAEAMETEES 



Model 








To 




CvO 






rad 


rad 




K (keV) 




erg cm"*^ 


erg cm~"^ 


Crab-like 


189 


2.4 X 10" 


-9 


2.1 X 10*^ 


(18) 


9.5 X 10^^ 


5.7 X 10^1 


Vela-liko 


70.4 


1.0 X 10" 


-10 


6.1 X 10^ 


(5.3) 


2.8 X 10^^ 


1.9 X 10^2 


Older star" 


10.5^ 


3.2 X 10" 


-13 


6.1 X 10^ 


(5.3) 


2.8 X 10^^ 


1.9 X 10^2 



"Spin-down age 5 X 10'-' years. 

^Scaled from Vela's spin rate with an n = 3 breaking index model. 



TABLE 2 

MODEL SIMULATION PAEAMETEES 



Rc (km) 


I (g cm-) 


p ( g cm ■^) 


(3 (rad s"^ K) 


A (m) 


h (km) 


(7 (m) 


vq (cm s ^) 


9.9 


8.6 X 10'''' 


1.5 X 10^" 


2.0 X 10^ 


200 


2.8 


20 


lO'^ 
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